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Abstract 

Updated stellar rates for the reaction ^^Mg(p,7)^^Al are calculated by using 
all available experimental information on ^^Al excitation energies. Proton 
and 7-ray partial widths for astrophysically important resonances are derived 
from shell model calculations. Correspondences of experimentally observed 
Al levels with shell model states are based on application of the isobaric mul- 
tiplet mass equation. Our new rates suggest that the ^^Mg(p,7)^^Al reaction 
influences the nucleosynthesis in the mass A>20 region during thermonuclear 
runaways on massive white dwarfs. 
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I. INTRODUCTION 



Explosive stellar burning of hydrogen in the mass A>20 range is characterized by a large 
number of proton capture reactions and /5-decays. The resulting network of nuclear processes 
is called the rp-process 0]. This process might be responsible for the energy production and 
nucleosynthesis in a variety of astrophysical sites with different temperature and density 
conditions. For example, in novae, typical peak temperatures range from Tg = 0.2 — 0.4 
0, with Tg the temperature in GK. For x-ray bursts and accreting black holes the rp- 
process could take place at very high temperatures in excess of Tg=l Stellar rates of 
several proton capture reactions relevant for the rp-process were estimated by Wallace and 
Woosley and by Wiescher et al. f^. Recently, some of these rates were updated with new 
experimental information and improved theoretical models 

At low stellar temperatures Tg<0.1 the isotope ^^Mg is synthesized in the NeNa-cycle. 
Under such conditions the /5-decay of ^^Mg (Ti/2=ll-3 s) and the subsequent ^^Na(p,a) 
reaction convert material back into ^°Ne, giving rise to cycling of material in the NeNa mass 
range. If the stellar temperature is sufficiently high the proton capture reaction on ^^Mg 
becomes faster than the competing /3-decay. In this case the reaction flow breaks out of 
the NeNa mass region and a whole range of heavier nuclei could be synthesized, depending 
on the temperature-density conditions and the duration of the astrophysical event. This 
scenario, for example, might be responsible for the synthesis of elements such as Si, S and 
Ar, which have been found to be overabundant in the ejecta of ONeMg novae Therefore, 
a quantitative estimate of the stellar reaction rate for ^^Mg(p,7)^^Al is important in order 
to model the nucleosynthesis in the mass A>20 range. At very high temperatures above 
Tg=l the ^^Mg(p,7)^^Al reaction is of minor importance, since the isotope ^^Mg is bypassed 
via the sequence 2iNa(p,7)22Mg(p, -ff^A\{p,-fy'^Si{f3+iyf'^Al 

The reaction rate for ^^Mg(p,7)^'^Al was previously estimated by Wallace and Woosley 

and their calculation was based on a single resonance only. Subsequently the reaction 
rate calculation was improved by Wiescher et al. |^ who considered three resonances and 
in addition a contribution from the direct capture process. The most recent estimate was 
pubhshed by Kubono et al. |]TD[. These authors, who investigated the level structure of ^^Al 
near the proton threshold by using the ^^Mg(^He,t)^^Al charge-exchange reaction, based 
the reaction rate estimate on their experimentally determined ^^Al excitation energies and 
spin-parity restrictions. 

We present a reanalysis of the ^^Mg(p,7)^^Al reaction rate for several reasons. First, a 
recent experimental study of the ^^Mg(^He,t)^^Al charge-exchange reaction was published by 
Greenfield et al. [|lT|. The excitation energies of proton threshold levels in ^^Al measured by 
the two groups differ by about 30-50 keV. This difference might change the resulting reaction 
rates appreciably. Second, the authors of Ref. |jlO| conclude that the analog assignments of 
the two lowest-lying proton threshold states in ^"^Al are still uncertain, resulting in large 
errors of the derived stellar reaction rates. Third, the proton and 7-ray partial widths of the 
resonances in question have never been measured. These quantities were crudely estimated 
in Ref. |]10[ by adopting 'typical' single-particle spectroscopic factors and 'average' 7-ray 
transition strenghts from Ref. 

In this work we use ^^Al excitation energies recommended by Ref. |]T2[ which are based 
on previously published experimental results. We present additional support for the analog 
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assignments of the proton threshold levels in ^''Al. Furthermore, we calculate the proton and 
7-ray partial widths of astrophysically important levels by using the nuclear shell model. In 
Sect. II we describe briefly the formalism for calculating stellar reaction rates. Experimental 
and theoretical nuclear input parameters are presented in Sect. III. Results and astrophysical 
implications are discussed in Sect. IV. A summary is given in Sect. V. 



II. CALCULATION OF STELLAR REACTION RATES 

The proton capture cross sections on sd-shell nuclei are predominantly determined by 
summing the contributions from isolated resonances corresponding to unbound compound 
nuclear states and from the nonresonant direct capture (DC) process. In the following we 
briefly describe the method of calculating the resonant and nonresonant (DC) contributions 
to the stellar reaction rates. 



A. Resonant Reaction Contributions 

For the reaction under consideration here the resonances are narrow and isolated. The 
resonant rate contribution can be calculated from resonance energies Ej and resonance 
strengths w-^i (both in units of MeV) [|13| 

A^A < (TV >,= 1.54 X 10^^(/ir9)"^/^^(cJ7)iexp {-11.605Ei/Tg) cm^^mole'^s"^ . (1) 

i 

The resonance strength for a (p,7) reaction is given by 

2 7-1-1 F r 

where J and jt are the spins of the resonance and the target nucleus, respectively, and the 
total width Ftot is the sum of the proton partial width Fp and the 7-ray partial width F^. 

The proton partial width Fp can be estimated from the single-particle spectroscopic 
factor S and the single-particle width Fgp of the resonance by using 

Fp = C25-Fsp , (3) 

where C is the isospin Clebsch-Gordan coefficient. Spectroscopic factors S are calculated in 
this work by using the nuclear shell model (Sect. III). Single-particle widths Fsp are obtained 
from resonant scattering phase shifts generated by an appropriate folding potential (see 
below). In this context, the partial width Fgp is defined as the energy interval over which 
the resonant phase shift varies from tt/4 to Sir/ 4. 

Gamma-ray partial widths for specific electromagnetic transitions are expressed in terms 
of reduced transition probabilities B( Jj ^/;L) which contain the nuclear structure informa- 
tion of the states involved in the transition |[T^. In the present work, the reduced transition 
probabilities are calculated in the framework of the nuclear shell model (Sect. III). The total 
7-ray width F^ of a particular resonance is given by the sum over partial 7-ray widths for 
transitions to all possible lower-lying nuclear states. 
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B. Nonresonant Reaction Contributions 



The nonresonant proton capture cross section is calculated by using the direct capture 
(DC) model described in |T^-[18|. The DC cross section af"^ for a particular transition 
is determined by the overlap of the scattering wave function in the entrance channel, the 
bound-state wave function in the exit channel and the electromagnetic multipole transition 
operator. Usually, only the dominant El transitions have to be taken into account. Wave 
functions are obtained by using a real folding potential given by [0,0] 



V{R) = A VpiR) = ^ J J Pa(ri)pA(r2) v^s (E, p„, pa, s) dridrg 



(4) 



Here A represents a potential strength parameter close to unity, and s = |R + r2 — ri|, with 
R the separation of the centers-of-mass of the projectile and the target nucleus. The mass 
density distributions pa and pA are either derived from measured charge distributions 
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or calculated by using structure models (e.g., Hartree-Fock calculations). For the effective 
nucleon-nucleon interaction Ves we used the DDM3Y parametrization |T^. The imaginary 
part of the potential has been neglected due to the small flux into other reaction channels. 

The total nonresonant cross section a^j- is determined by summing contributions of direct 
capture transitions to all bound states with single-particle spectroscopic factors 



DC 



(5) 



The astrophysical S-factor of a charged-particle induced reaction is defined by |]T3 

S{E) = Eexp {27vr])a{E), 



(6) 



with 1] denoting the Sommerfeld parameter. If the S-factor depends only weakly on the 
bombarding energy the nonresonant reaction rate as a function of temperature Tg can be 
expressed as 



A^A < (TV >j 



7.833 X 10' 



X exp 



-4.249 




5(Eo)[MeVbarn] 



cm mole s 



i„-i 



(7) 



with Zi and Z2 the charges of the projectile and target, respectively, and A the reduced 
mass (in amu). The quantity Eq denotes the position of the Gamow peak corresponding to 
the effective bombarding energy range of stellar burning. 



III. EXPERIMENTAL AND THEORETICAL INPUT PARAMETERS 

In this section we present a discussion of excitation energies, spectroscopic factors and 
7-ray partial widths which enter in the calculation of stellar reaction rates. 

Experimental excitation energies in ^^Al below E^=3 MeV have been compiled by Endt 
!l[|. The energies listed in Table 24.23 of Ref. |^ are based on ^^Mg(^He,t) and (p,n) 
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reaction studies performed prior to 1990. Three recent charge-exchange reaction studies 
1^,11;^ also report ^"^Al excitation energies. In the present work we use the Ea.(^^Al) 



values compiled and evaluated by Ref. |T2[ which are based on experimental information. 



The values measured in the (p,n) study of Kiang et al. p3| have been disregarded because 
of the superior energy resolution and counting statistics of the ('^He,t) reaction studies. Our 
adopted excitation energies (Tables I and II) differ from the results of Kubono et al. on 
average by about 20 keV. 

Spectroscopic factors and reduced 7-ray transition strengths for the levels of astrophys- 
ical interest are calculated by using the nuclear shell model. This procedure requires the 
identification of experimentally observed ^^Al levels with calculated shell model states. For 
bound ^^Al states there is a one-to-one correspondence. However, the spins and parities of 
experimentally observed unbound states are not known uniquely, resulting in ambiguities 
for the shell model assignments. Level assignments based on a comparison of experimental 
and shell model excitation energies alone are not useful either, mainly due to the difficulty 
of the shell model in producing accurate E^. values. In this work we have used a method 
described in Ref. to which the reader is referred for details. In brief, experimental excita- 
tion energies of ^^Na and ^^Mg states, for which the spins and parities are well-known, are 
used together with the isobaric multiplet mass equation (IMME) in order to calculate 
excitation energies of ^^Al analog states: 

i?.rAl) = 2E.rMg)-E.rNa) + 2[c-c(g.s.)] . (8) 

The coefficient c is a measure for the isotensor Coulomb energy of a specific isobaric triplet 
and is estimated in this work by using the nuclear shell model (see below). Correspondences 
between experimentally observed ^^Al levels and shell model states are found by a) minimiz- 
ing the difference between measured excitation energies and E^. values calculated from Eq. 
(^), and b) matching experimentally determined spin-parity restrictions with shell model 
quantum numbers. It can be seen from our results listed in Table | that the experimentally 
observed ^^Al states at Ea;=2349 and 2534 keV most likely correspond to the 83" and 



shell model states, respectively, in agreement with the tentative assignments of Ref. pO 
However, the experimental states at Ea;=2810 and 2900 keV most likely correspond to the 
2j and 3^ shell model states, respectively, in contradiction with the results of Kubono et 



al. 10 



Shell model calculations have been performed by using the code OXBASH The 
isospin-nonconserving (INC) interaction of Ormand and Brown [^] is employed for the 
calculation of wave functions and excitation energies of T=l triplet states in the mass 
A=24 system. Coefficients c in Eq. (8) are estimated from theoretical excitation energies. 
Shell model wave functions are used for calculating single-particle spectroscopic factors and 
reduced transition probabilities for Ml- and E2- 7-ray decays. Resulting values of the 
resonance parameters are listed in Table ||. It can be seen that the 7-ray partial width is 
much smaller than the proton partial width for all resonances considered. Therefore, the 
resonance strength depends mainly on the value of F^. Table |Tl| displays the calculated 
F^ values for the ^^Al unbound states of main interest in this work together with theoretical 
and experimental F^ values of the corresponding ^^Na mirror states. It can be seen that in the 
case of ^^Na the shell model 7-ray widths are in excellent agreement with the experimental 
values deduced from lifetime measurements. 
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In Table II our resonance strengths are compared to previous results [|T0|- For the first 
resonance the strengths are very similar, since the experimentally measured lifetime of the 
^^Na mirror state was used in Ref . . For the second resonance our strength a;7 is consid- 
erably smaller than the value of Kubono et al. [l^ who adopted an 'average' 7-ray transition 
strength from Ref. For the third and fourth resonance the spin-parity assignments are 
interchanged compared to Ref. ||lOl (see above) which explains the discrepancy of the 077 
values. It should be noted that the latter two resonances have a negligible influence on the 
stellar reaction rates (Sect. IV). 

The parameters for the direct capture (DC) contribution to the stellar reaction rates are 
presented in Table |IV|. All transitions considered are displayed together with our calculated 
shell model spectroscopic factors. With this information, the total astrophysical S-factor 
for the direct capture process into all bound states has been determined (Sect. lib). For 
bombarding energies below 1 MeV the S-factor can be expressed as 



S{E) = 22.5 - 1.1 X IQ-^E + 6.9 x lO'^E^ keV ■ b. 



(9) 



Our derived direct capture S-factor is about 10 % smaller than the results of Wiescher et 
al. |H] which were also adopted by Ref. 110]. 



IV. DISCUSSION 

The recommended stellar rates of the ^^Mg(p,7)^^Al reaction can be parametrized for 
temperatures below T9=2 by the expression 

Na <av>=J2 A/T932 exp {-Bi/T^) + C/T923 exp (-LI/TQIS) cm^mole'^s"^ , (10) 

i 

where, for example, T932 stands for Tg^^. The first and second term in Eq. (|T0|) represent 
the contributions of all narrow resonances and the direct capture process, respectively. The 
parameters Aj, Bi, C and D are listed in Table ^ 

The various contributions to the total reaction rate are displayed in Fig. |l|. The direct 
capture process determines the stellar rates at low temperatures of Tg < 0.2. The E/j=478 
keV resonance dominates the reaction rates in the range Tg=0. 2-1.0. The E^=663 keV res- 
onance is of importance at high temperatures above Tg=l only. The resonances at Eij=939 
and 1029 keV are negligible over the whole temperature range. Our results are compared 
in Fig. with previous work At stellar temperatures above Tg=l the reaction rates 



of the present work are smaller than the results of Refs. ||5|JT0|1 by about 70%. In the tem- 



perature range Tg=0.2-0.5 important for hydrogen burning in novae (see below) the present 
reaction rates deviate up to a factor of 3 from the values given in Ref. [1^, and up to a 
factor of 2 from the results of Ref. P|. The reaction rates for ^^Mg+p are therefore now 
based on more consistent experimental and theoretical input parameters. 

Figure 3 presents temperature and density conditions for which the proton capture re- 
action on ^^Mg and the ^^Mg /3-decay are of equal strength. The solid line is calculated by 
assuming a hydrogen mass fraction of X//=0.365 [0. Recent results of hydrodynamic studies 
of ONeMg novae are also shown in Fig. 3. The full circles represent temperature and 
density conditions at the peak of the thermonuclear runaway for accretion onto white dwarfs 
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of different inital masses (I.OOMq, 1.25M0 and l.SSM©). Our results indicate that for white 
dwarfs of masses <1.25M0 the proton-capture reaction on ^^Mg is slower than the competing 
/3-decay and, therefore, is of minor importance for the resulting nucleosynthesis. However, 
for accretion onto very massive white dwarfs (I.SSMq model of Ref. the ^^Mg(p,7)^'^Al 
reaction dominates over the ^^Mg /5-decay and will influence the nucleosynthesis in the mass 
A>20 range. Stellar network calculations are underway in order to investigate quantitatively 
the implications of our new ^^Mg+p reaction rates and their corresponding uncertainties. 



The results will be published in a forthcoming paper 26 



V. SUMMARY AND CONCLUSIONS 



Improved estimates of reaction rates for ^^Mg(p,7)^^Al are presented in this work. Res- 
onance energies are derived from all available experimental information on ^^Al excitation 
energies. Proton and 7-ray partial widths of astrophysically important resonances are esti- 
mated from single-particle spectroscopic factors and reduced 7-ray transition probabilities, 
respectively, derived from shell model calculations. Correspondences of experimentally ob- 
served ^^Al levels with shell model states are based on application of the isobaric multiplet 
mass equation. In the temperature range T9=0.2-0.5 important for the nucleosynthesis 
in novae the present reaction rates deviate up to a factor of three from previous results. 
Our new stellar rates suggest that the ^^Mg(p,7)^'^Al reaction will influence the nucleosyn- 
thesis in the mass A>20 region during thermonuclear runaways on massive white dwarfs. 
Quantitative predictions have to await the results of large-scale stellar network calculations. 
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TABLES 



TABLE I. Comparison of excitation energies (in MeV). 





94 A ^n 


9Atvt 

^*Na 




24 Al (IMME) ^ 


24 Al (OXBASH) ^ 


Ex 


J"" 


Ex 


J"" 


Ex 






J"" 


0.0 


4+ 


0.0 


4+ 


0.0 


4+ 


0.0 


4+ 


0.426 


1+ 


0.472 


1+ 


0.441 


1+ 


0.448 


1+ 


0.510 


2+ 


0.563 


2+ 


0.450 


2+ 


0.580 


2+ 


1.107 


(1-3)+ 


1.347 


1+ 


1.073 


1+ 


1.100 


1+ 


1.130 


fl -3)+ 


1.341 


2+ 


1.069 


2+ 


1.126 


2+ 


1.275 


3+ 


1.345 


3(+) 


1.283 


3+ 


1.374 


3+ 


1.559 


(5+) 


1.512 


5+(3+) 


1.552 


5+ 


1.554 


5+ 


1.559 


(2+) 


1.846 


2+ 


1.541 


2+ 


1.590 


2+ 


1.634 


3+ 


1.886 


3+ 


1.715 


3+ 


1.729 


3+ 


2.349 


(3+) 


2.514 


3+ 


2.305 


3+ 


2.176 


3+ 


2.534 


(4,5)+ 


2.563 


4+(2+) 


2.476 


4+ 


2.541 


4+ 


2.810 


2+ 


2.978 


2+(3+) 


2.803 


2+ 


2.837 


2+ 


2.900 


(1-3)+ 


2.904 


3+ 


2.737 


3+ 


2.629 


3+ 



Experimental values evaluated and compiled in Ref. fl^ . 

Calculated by using Eq. (8). 

Calculated by using the nuclear shell model. 



TABLE IL Resonance parameters for the reaction Mg(p,7) Al 



E^ (MeV) " 




(MeV) ^ 


Lp (meV) 


(meV) 


LOj (meV) 


ioj (meV) 


2.349±0.020 




0.478 


185 


33 


25 


27 


2.534±0.013 




0.663 


2.5 X 10^ 


53 


58 


130 


2.810±0.020 


2+ 


0.939 


9.5 X 10^ 


83 


52 


11 


2.900±0.020 


3+ 


1.029 


3.4 X 104 


14 


12 


16 



12 



" Experimental values adopted from Ref, 

Calculated from column 1 and Qp^=1871±4 keV p5 |. 
= From Ref. 0. 



9 



TABLE III. Gamma-ray partial widths (in meV) of ^^Al and ^^Na levels. 





(24 Al) a 


Tf^ (24Na) « 


pexp (24]^g_) 6 




33 


43 


46 ± 14 




53 


51 


> 27 




83 


93 


> 27 


3| 


14 


15 


13 ± 2 



" Shell model values for the ^^Al - ^^Na mirror pair. 
Determined from measured lifetimes of ^^Na levels. 



TABLE IV. Shell model spectroscopic factors of Al bound states. 



Ea.(MeV) 




p ld3/2 


c^s 

P ld5/2 


P 2si/2 


0.000 


4+ 




0.39 




0.426 


1+ 


0.0010 


0.69 


0.0031 


0.510 


2+ 


0.0000 


0.29 


0.052 


1.107 


1+ 


0.098 


0.016 


0.36 


1.130 


2+ 


0.084 


0.0001 


0.30 


1.275 


3+ 


0.024 


0.0057 




1.559 


5+ 








1.559 


2+ 


0.11 


0.088 


0.19 


1.634 


3+ 


0.0009 


0.19 





TABLE V. Recommended parameters for the ^^Mg(p,7)^^Al reaction rate." 



Ai Bi C D 

4.02 X 10^ 5^56 3.72 x 10^ 21.95 

9.59 X 10^ 7.71 

8.52 X 10^ 10.91 

2.0 X 10^ 11.95 



" The total stellar reaction rate is given by Eq. (10). 
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FIGURES 




0.1 1 

stellar temperature T (GK) 

FIG. 1. Total stellar rate (solid line) and individual contributions (dashed lines) for the reaction 
23Mg(p,7)24Al. 
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I I I 



T 1 — I I I 



Present/Wiescher 




Present/Kubono 



J 1__L 



J I I I I I 



0.1 

0.1 1 
stellar temperature T (GK) 

FIG. 2. Ratio of the present reaction rate to previous results of Wiescher et al. and Kubono 
et al. [10|. The reaction rates are based on measured ^^Al excitation energies (see text) which are 
most likely Gaussian distributed. 
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10" - 



10^ - 



^ 10 h 
c 
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-2 



— \ 


1 1 1 1 1 1 1 

2^Mg(p,7)2'^AI - 


— 


\ 1.35 


1.00 


\ • ~ 


• 


1.2S\ 


1 





0.1 1 
Temperature T (GK) 

FIG. 3. Temperature-density boundary at which the proton capture reaction on ^'^Mg and the 
^^Mg /?-decay are of equal strength, assuming a hydrogen mass fraction of Xj7=0.365 Q. Peak 
temperature and density conditions achieved in the nova models of Ref. [Q] are indicated by full 
circles (see text). 
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